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Ram pressure stripping of disc galaxies orbiting in clusters. 
I. Mass and radius of the remaining gas disc 
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ABSTRACT 

We present the first 3D hydrodynamical simulations of ram pressure stripping of a 
disc galaxy orbiting in a galaxy cluster. Along the orbit, the ram pressure that this 
galaxy experiences varies with time. In this paper, we focus on the evolution of the 
radius and mass of th e remaining ga s disc a nd compare it with the classical analytical 
estimate proposed by iGunn fc Gottl (|1972[ ). We find that this simple estimate works 
well in predicting the evolution of the radius of the remaining gas disc. Only if the ram 
pressure increases faster than the stripping timescale, the disc radius remains larger 
than predicted. However, orbits with such short ram pressure peaks are unlikely to 
occur in other than compact clusters. Unlike the radius evolution, the mass loss history 
for the galaxy is not accurately described by the analytical estimate. Generally, in the 
simulations the galaxy loses its gas more slowly than predicted. 

Key words: galaxies: spiral - galaxies: evolution - galaxies: clusters - intergalactic 
medium 



1 INTRODUCTION 

Galaxies populate different environments, ranging from iso- 
lated field regions to dense galaxy clusters. Depending on 
the environment, the properties of galaxies, especially disc 
galaxies, change: In denser regions, disc galaxies tend to 
contain less neutral gas, show a weaker star formation ac- 
tivity and_^edder colou r s than galaxies in sparse regions 
(e.g. Ivan Gorkoml 120041 : iGoto et al.1 120031 and references 
therein). 

Ram pressure stripping (RPS), i.e. the removal of a 
galaxy's gas disc due to its motion through the intraclus- 
ter medium (I CM) , is one of the like ly candidates to explain 
these features. IGunn fc Gottl (jl972h (GG72) proposed that 
for galaxies moving face-on through the ICM the success or 
failure of RPS can be predicted by comparing the ram pres- 
sure with the galactic gravitational restoring force per unit 
area. Hydrodynamical simulations of RPS presented so far 
(jAbadi et al.1 ll 9991. iQuilis et al.ll2000l ISchulz fc Struck 2001 



(SSO U jMarcolini et al.ll2003l fMBD03).lRoediger fc Hensler 
120051 fRH05). iRoediger fc Briiggenl liood fRB06l) suggest 
that this analytical estimate does indeed do a fair job as long 
as the galaxies are not moving close to edge-on. However, in 
all these simulations the model galaxy was not exposed to a 
varying ICM wind, as one expects for cluster galaxies, but 
to a constant one. A few 2D simulations of spherical galaxies 



on radial orbits in clusters ([Lea fc Youngl 19761 , iTakeda et al.l 
1984, lAcreman et aH [20Q3) exist , where the ear l y ones were 
run with low resolution. iToniazzo fc Schindlerl (|200lh pre- 
sented 3D hydro-simulations for a spherical galaxy on a clus- 
ter orbit. Vollmer et al. (e.g. 1999, 2000, 2001, 2001a, 2003) 
have simulated the cluster passage of disc galaxies with a 
sticky-particle code, where the ram pressure is modelled as 
an additional force on the disc gas particles. Their work 
shows the importance of the modelling of realistic orbits, 
e.g. to determine the stripping history of individual galax- 
ies. However, the sticky-particle code cannot model hydro- 
dynamical effects such as instabilities that have been shown 
to play a role by the hydro-simulations. Additionally, the 
work of SS01, RH05 and RB06 has shown that the gas re- 
moval from the galactic potential does not occur instanta- 
neously, but that it takes some time until the gas is acceler- 
ated enough to become unbound from the galaxy's potential. 
Thus, in cases of short ram pressure peaks, this time delay 
may play an important role. 
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Here we present the first 3D hydrodynamical simula- 
tions of RPS a of disc galaxy orbiting in galaxy cluster. We 
focus on the question, in how far a time-dependent version 
of the analytical estimate based on the GG72 criterion can 
describe the evolution of the remaining gas disc. 
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2 METHOD 



2.2 Model galaxy 



We model the flight of a disc galaxy through a galaxy cluster. 
The galaxy starts at a position ~ 1 to 2 Mpc (depending on 
orbit) from the cluster centre with a given initial velocity. 
We use analytical potentials for the galaxy and the clus- 
ter, as this reduces computational costs significantly. Given 
the high velocities of cluster galaxies, the tidal effect on the 
galaxy is expected to be small. The work of Moore et al. 
(1996, 1998, 1999) demonstrated that only harassment, i.e. 
the cumulative effect of frequent close high velocity encoun- 
ters between cluster galaxies and the overall tidal field of 
the cluster affects cluster galaxies seriously. Thus, we ex- 
pect that our treatment yields reasonable results. 

The orbit of the galaxy is determined by integrating the 
motion of a point mass through the gravitational potential 
of the cluster. In the course of the simulation, the position 
of the galaxy potential is shifted along this orbit. 



2.1 Code 

The simulations w ere performed with the FLASH code 
([Frvxell et al.ll2000h , a multidimensional adaptive mesh re- 
finement hydrodynamics code. It solves the Riemann prob- 
lem on a Cartesian grid using the Piecewise-Parabolic 
Method (PPM). The simulations presented here are per- 
formed in 3D. The gas obeys a polytropic equation of state 
with an adiabatic index of 7 = 5/3. The size of the sim- 
ulation box is chosen such that the galaxy's orbit during 
the simulation time (3 Gyr) fits into the grid. Depending 
on the orbit, the size of the simulation box ranges between 
(2 Mpc) 3 and 2x5x2 Mpc 3 . For an example see Fig. Q] All 
boundaries are reflecting. 

The coarsest refinement level has a resolution of Ax = 
62.5 kpc. For most runs, we use 8 levels of refinement, i.e. the 
best resolution is Ax = 0.5 kpc. In addition to the standard 
density and pressure gradient criteria, our user-defined re- 
finement criteria enforce maximal refinement on the galactic 
gas disc and enforce stepwise de-refinement with increasing 
distance to the galaxy in the following fashion: 

• We always refine a sphere of radius 50 kpc around the 
galaxy centre to at least Ax = 3.9 kpc (5 refinement levels). 

• We always use the maximum refinement level for the 
galactic disc region, which is defined as the cylinder with 
radius 27 kpc and thickness 6 kpc around galactic centre, 
where the density larger than 10 -27 gem -3 . 

• We limit the refinement outside 50 kpc around galaxy 
centre to Ax = 2 kpc (6 refinement levels). 

• We limit the refinement outside 150 kpc around galaxy 
centre to Ax = 15.6 kpc (3 refinement levels). 

• Additionally, we always refine a sphere of radius 100 kpc 
around cluster centre to Ax = 7.8 kpc (4 refinement levels). 

Our user-defined de-refinement rules overwrite the standard 
pressure and density gradient criteria, but our user-defined 
refinement criteria overwrite the de-refinement criteria. The 
FLASH code also ensures that neighbouring blocks differ by 
maximal one refinement level. 

A detailed discussion of the influence of the resolution 
on our results is given in Sect. [Al 



The galaxy model is the same as in RB06, i.e. a massive 
spiral with a flat rotation curve at 200 km s -1 . It consists 
of a dark matter halo, a stellar bulge, a stellar disc and a 
gaseous disc. All non-gaseous components just provide the 
galaxy's potential and are not evolved during the simulation. 
For a description of the individual components and a list of 
parameters please refer to RB06. 

Initially, the galaxy is set in pressure equilibrium with 
the surrounding ICM. For the case of a homogeneous ICM 
density and pressure, this task is straightforward. In our 
case, however, there is the difficulty that both, pressure and 
density of surrounding ICM, vary with distance from the 
cluster centre. Consequently, the ICM pressure at the "sur- 
face" of the galaxy is not the same everywhere. Therefore, 
we use the following procedure for setting the galaxy in pres- 
sure equilibrium with the surrounding ICM: 

(i) We chose values for the ICM density and pressure 
somewhat lower than the values at the galaxy's initial posi- 
tion. 

(ii) We calculate the density, pressure and rotation ve- 
locity in the galaxy's gas disc such that the gas disc is in 
pressure equilibrium with a homogeneous ambient gas of 
density and pressure chosen in (i). 

(iii) We now need to ensure a continuous transition from 
the galactic disc to the cluster ICM. Therefore, we compare 
the pressure calculated in (ii) with the ICM pressure in each 
grid cell. Inside the gas disc the pressure calculated in (ii) is 
larger than the ICM pressure. However, at the transition to 
the ambient gas used in (ii), the pressure set in (ii) will be 
lower than the ICM pressure. In such grid cells, we substitute 
the values from (ii) with ICM values. This ensures that the 
pressure distribution is continuous everywhere. However, the 
density distribution shows a jump between the ICM and the 
ISM. Rotation velocity plus initial orbital velocity is only 
set for ISM cells. 

Figure [2] displays radial profiles for the density, projected 
ISM surface density, pressure, and total velocity in the galac- 
tic plane for two runs. The galaxy's initial position in these 
two runs differs (see Table [2}. Therefore, the gas discs have 
to be in pressure equilibrium with different external pres- 
sures and thus differ slightly in radius. Moreover, the radius 
of each gas disc varies slightly with azimuthal angle, because 
the pressure in the ICM surrounding the galaxy is not homo- 
geneous. However, at least in the inner 10 kpc, the profiles 
of different runs are indistinguishable. 

2.3 Cluster model 

The cluster is set in hydrostatic equilibrium, where the den- 
sity of the ICM follows a /^-profile, 

-3/2/3 



Picm{t) — piCMO 



1 + 



Ric 



(1) 



while the ICM temperature is constant. Here we present 
simulations for three clusters. The parameters are listed in 
Table [1] Figure [3] shows the density and pressure profiles 
of clusters CI and C3. The only difference between clusters 
CI and C2 is a factor of 2 in the density. These two clusters 
are compact, the ICM density is strongly peaked. Thus, they 
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Figure 1. Colour-coded gas density in the orbital plane for run C2-SM-SLW-MFM (here with best resolution 0.25 kpc). The lhs picture 
shows a slice through the complete simulation box, the rhs a blow-up of the region around the galaxy. The black line marks the galaxy's 
orbit. The white rectangles show the adaptive grid: Each rectangle marks the boundaries of one block, each block consists of 8 3 grid 
cells. 



Table 1. ICM parameters: Core radius, Ricm, beta-parameter, 
(3, central ICM density, Picmo, and ICM temperature, Ticm> f° r 
clusters CI, C2 and C3. 



#iCM/kpc 


P 


PlCMo/(gcm 3 ) 


Ticm/ K 


CI 50 


0.5 


2 • 10- 26 


4.7- 10 7 


C2 




1Q -26 




C3 386 


0.705 


6.07- 10- 27 


9.5 • 10 7 



Mohr et alJll999h , which 



are similar to the Virgo cluster ([Matsumoto et al1 l2000). but 
not as extreme as Virgo (see Fig. [3]) . Cluster C3 r esembles 
the Coma cluster (parameters see 
is very extended. 



2.4 Galaxy orbit 

The orbit of the galaxy determines its ram pressure his- 
tory. According to the simulations of RB06, for our model 
galaxy, complete stripping is expected for ram pressures of 
about 10 -10 erg cm -3 . We will refer to such ram pressures 
as "high". "Medium" ram pressures (~ 10 -11 erg cm -3 ) are 
still expected to strip the gas disc significantly. We aimed 
at constructing orbits with medium to high ram pressure 
peaks. Moreover, we concentrate on orbits of galaxies that 
could be regarded as falling into the cluster for the first 
time, i.e. they start from a sufficiently large distance from 
the cluster centre. Additionally, the galaxies are bound to 
the cluster. 

Figure |4] shows several representative orbits in clusters 
CI and C3. Already at this stage some important issues can 
be learned: 

• Naturally, orbits with large impact parameters yield 
ram pressure peaks of long duration. Short ram pressure 
peaks are typical for galaxies passing near the cluster cen- 
tre. However, even on orbits with small impact parameters, 



the ram pressure peaks in the extended cluster are about 
twice as wide as on comparable orbits in the compact clus- 
ter. 

• Orbits with low ram pressure peaks are practically im- 
possible in both, the extended and the compact cluster. Even 
on orbits with large impact parameter, where the ICM den- 
sity along the orbit is small, the galaxy's velocity is high 
enough to yield at least medium ram pressures. If the galaxy 
passes near the cluster centre, typically both, ICM density 
and galaxy velocity are high. 

• In the extended cluster, even orbits with medium ram 
pressure peaks are difficult to construct. Almost all galaxies 
should be stripped severely or completely. 

• Comparing the pressure profile along the galactic disc as 
shown in Fig. [2] with the radial pressure profiles of the clus- 
ters (Fig. [3} gives an interesting insight into the importance 
of the cluster environment for the galaxy: The ICM-pressure 
in the Coma-like cluster drops below 10 -12 erg cm -3 only for 
radii larger than ~ 2 Mpc. Inside ^ 1 Mpc the ICM pressure 
is above 10 -11 erg cm -3 . These values are comparable to the 
pressure in the galactic plane at disc radii of approx. 15 and 
8 kpc. Thus, even if a galaxy could reach such a position in 
the cluster without suffering from any other processes, its 
gas disc should be affected at least by the external pressure. 

• During most of the orbit, the ram pressure is even larger 
than the thermal ICM pressure (bottom panel of Fig. [SJ. 



Here we present simulations for four orbits that are sum- 
marised in Fig. [5] The labels of the runs represent the cluster 
(C+number), small or large impact parameter (SM or LG), 
fast or slow galaxy (FST or SLW), inclination (F for near 
face-on, M for medium, E for near edge-on; FE for first near 
face-on but then near edge-on, etc.). In all cases, the galaxy 
orbits in the y-z-plsme. Table [2] lists the simulation runs. 
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Figure 2. Radial profiles of the density, p, pressure, p, and total 
velocity, vtot, in the galactic plane for the initial model for two 
different runs (see legend). Also the radial profile for the projected 
ISM surface density, S, is shown. One dot for each cell in the 
galactic plane is plotted. In the bottom panel, the total velocity in 
the galaxy's rest frame is shown. Velocities inside the disc radius 
(as set by e.g. the density profile) are due to the galaxy's rotation, 
while velocities outside the disk radius are due to the galaxy's 
motion through the ICM. 
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Figure 3. Density and pressure profiles for model clusters CI 
and C3. The profiles for cluster C2 are a factor of 2 lower than 
the ones of cluster CI. For comparison, the profiles for the Virgo 
cluster are shown. 



Table 2. List of runs. Run labels code the cluster (C+number), 
small or large impact parameter (SM or LG), fast or slow galaxy 
(FST or SLW), inclination (F for near face-on, M for medium, E 
for near edge-on; FE for first near face-on but then near edge-on, 
etc.). Initial galaxy coordinates £ ga io and z ga io are always zero, 
2/gal0 is given in the second column. The third column lists the 
initial galaxy velocity. The inclination listed in the fourth column 
is the angle between the galaxy's rotation axis and the y-axis. 



label 


2/gal0 


initial velocity 


inclination 




in kpc 


in km s _1 


in ° 


Cl-SM-FST-MF 


1500 


(0, -600, 100) 


-30 


Cl-SM-FST-E 






80 


Cl-LG-FST-MF 


1300 


(0, -600, 500) 


-30 


Cl-LG-FST-EF 






-60 


Cl-LG-FST-FE 






45 


C2-SM-SLW-FME 


800 


(0,0,200) 


30 


C2-SM-SLW-EMF 






-60 


C2-SM-SLW-MFM 






-20 


C3-SM-FST-MF 


2300 


(0,-800,200) 


30 



3 RESULTS 

Previous simulations where the galaxy was exposed to a 
constant ICM wind found that the ICM-ISM interaction is 
a multi-stage process (SS01, MBD03, RH05, RB06). The 
stages reflect the two gas removal mechanisms that work on 
different timescales. 

One process is ram pressure "pushing", the immediate 
consequence of the ram pressure dislocating the gas disc at 
radii where the galaxy's gravity is not strong enough. This 
process results in the instantaneous stripping phase where 
the gas is pushed out of its original position on a timescale 
of a few 10 Myr. However, the gas does not become unbound 
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Figure 4. Representative galaxy orbits in clusters CI and C3: 
The top panel shows the orbits, which are in the y-z-plane. The 
next panels show the temporal evolution of ram pressure, p ram , 
galaxy velocity, v ga i, ICM density along orbit, picM? distance 
to cluster centre, d, and ratio between ram pressure and local 
ICM pressure, Pram /PiCM- The clusters (CI or C3), typical galaxy 
velocity (FaST or SLoW) and impact parameter (SMall or LarGe) 
are coded by line colour and thickness, see legend. 
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Figure 5. Summary of galaxy orbits for our simulations. The 
top panel shows the orbits, which are in the y-z-plane. Crosses 
mark the position of the galaxy in intervals of 250 Myr. The next 
panels show the temporal evolution of ram pressure, p ram , galaxy 
velocity, v ga \, ICM density along orbit, picm, distance to cluster 
centre, d, and ratio between ram pressure and local ICM pressure, 
Pram/piCM- For an explanation of the labels see text. 
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from the galactic potential immediately but with a certain 
delay, because it takes some time until the gas is accelerated 
to the escape velocity This time delay grows with decreasing 
ram pressure as the ram pressure is responsible for acceler- 
ating stripped gas packages. 

The second gas removal process is con tinuous strip- 
ping (or turbulen t /viscous stripping, see e.g. Nulsen 1982; 
iQuilis et al. 2000), which is caused by the Kelvin- Helmholtz 
instability induced by the ICM wind flowing over the sur- 
face of the gas disc. Continuous stripping leads to a slow 
but continuous gas loss at a rate of ~ IMq yr _1 , it works 
on a timescale of Gyrs. The gas disc radius is not changed 
significantly by continuous stripping. 



3.1 Analytical estimates 

The usual way to estimate the amount of gas loss due to 
ram pressure pushing for galaxies moving face-on follows the 
suggestion of GG72. Here, one compares the gravitational 
restoring force per unit area, 



/grav(i?) — Hgas(i?) ^ W > 



and the ram pressure, 



Pram = PlCM^galj 



(2) 
(3) 



for each radius of the galaxy (see e.g. RB06), where <I> is the 
gravitational potential of the galaxy, S gas (R) the gas surface 
density, /?icm the ICM density and v ga \ the galaxy's velocity. 
At radii where the restoring force is larger, the gas can be 
retained, at radii where the ram pressure is larger, the gas 
will be stripped. The transition radius is called the stripping 
radius. In order to apply this estimate to our simulations, 
we compare the gravitational restoring force to the current 
ram pressure at every time. This estimate should work best 
for galaxies moving face-on. 

(Nearly) pure continuous stripping affects galaxies mov- 
ing edge-on. All galaxies with other inclinations experience 
a mixture of ram pressure pushing and continuous stripping. 
For a spherical gas cloud, the mass loss rate d ue to continu- 
ous stripping was estimated by iNulsenl (|1982). For a sphere 
of radius, R, moving though the ICM of density, picm, with 
velocity, v s& i, the appropriate mass loss rate is 7i\R 2 /9icM^ g ai- 
As the surface area of a sphere is twice as large as the surface 
area of a disc (counting upper and lower side), we use 



Mkh — 0.5nR picM^gai 



(4) 



for our case of a disc galaxy. Both, ICM density, picm, and 
galaxy velocity, v ga \ are time-dependent. As disc radius, R, 
we use the current analytical stripping radius as estimated 
by the GG72 criterion, which is also time-dependent. This 
choice may seem mismatched at first, but it is motivated by 
the fact that the GG72 criterion gives good results for the 
stripping radius for most inclinations. 

3.2 Snapshots 

Figure [6] shows snapshots for run C2-SM-SLW-MFM (high 
resolution). At larger distances to the cluster centre, the 
galaxy moves subsonically, in the inner cluster it moves su- 
personically as indicated by the bow shock. Before pericentre 
passage, the gas disc size decreases, whereas it stays rather 



constant after pericentre passage. Snapshots for other runs 
are similar. 



3.3 Comparison with numerical result 

Figure [3 compares analytical estimates and numerical re- 
sults for the stripping radius and the remaining gas mass. On 
the way towards the cluster centre, the distinction between 
the instantaneous and continuous stripping phases vanishes 
as both processes overlap. In earlier simulations with con- 
stant ICM wind (RH05, RB06, but also MBD03) a short 
period of gas backfall during the intermediate phase was 
observed. The gas disc mass as a function of time showed a 
characteristic local maximum after the instantaneous strip- 
ping phase. Now, in the simulations with varying ICM wind, 
we observe a repeated backfall of stripped gas prior and 
during pericentre passage. The backfall appears as repeated 
local maxima or plateaus in the gas disc mass as a func- 
tion of time. Intere stingly, the sticky-particle simulations of 
IVollmer et al.l ([200 lh where the galaxies were also exposed to 
a varying ram pressure found gas backfall of up to 5 x 10 8 Mq 
after pericentre passage, primarily for edge-on galaxies. This 
behaviour cannot be observed in our simulations. The back- 
fall cycles in our simulations are caused by the hydrodynam- 
ical velocity field in the galaxy's wake. On the way out of 
the cluster, the instantaneous stripping is not active any- 
more because all gas from outer radii is already lost. Here 
continuous stripping dominates. 



3. 3. 1 Ram pressure pushing 

The result of the analytical estimate that only takes ram 
pressure pushing into account is shown as black lines in 
Fig. [7] The analytical and numerical result for the strip- 
ping radius agree remarkably well. Only during the pericen- 
tre passage in runs C2-SM-SLW-. . . , the simulated galaxies 
can retain a larger gas disc than predicted. Here, the ram 
pressure, pramj increases faster than the stripping timescale: 
During the ~ 200 Myr before the ram pressure maximum, 
the stripping radius is predicted to decrease from 10 kpc to 
4kpc, which corresponds to an average radius decrease of 
1 kpc per 35 Myr. A rough estimate of the acceleration of 
gas with a surface density, E, is p r am/E. Thus, the stripping 
timescale, i.e the time needed to move the gas package a 
distance, s, is 



r — 




£ 10 -11 ergcm -3 
5 kpc 10 _3 g cm -2 



Pram 



(5) 



where we have inserted typical values for this case. How- 
ever, this estimate is a lower limit, as it neglects the effect 
of the galaxy's restoring force, which lessens the effective ac- 
celeration. Although this timescale is rather short, in runs 
C2-SM-SLW-. . .the ram pressure changes on an even faster 
timescale during pericentre passage. 

For all other runs, deviations of the gas disc radius from 
the analytical result can be attributed to different inclina- 
tions and they agree with the results of constant ICM wind 
simulations of RB06. As expected from the simulations of 
RB06, with a increasing deviation from the face-on inclina- 
tion, the remaining gas disc becomes increasingly asymmet- 
rical. This behaviour translates into an increasing distance 
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Figure 6. Snapshots for run C2-SM-SLW-MFM (high resolution): colour-coded gas density in the orbital plane. The black line marks 
the galaxy orbit. 
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Figure 7. Comparison between analytical and numerical stripping radius (top panels) and gas mass (middle panels). The bottom panels 
display the evolution of the true inclination, i.e. the angle between the galaxy's rotation axis and the direction of motion. Each subfigure 
is for one orbit (see title), different inclinations are colour-coded (see legend). For the numerical stripping radius, the mean value (thick 
lines) as well as maximum and minimum radius (thin lines) as described in RB06 are shown. For the gas mass, we show the bound gas 
mass (thick lines) and the mass of gas in the disc region (thin lines of matching colour, see also RB06). Analytical estimates according 
to the Gunn&Gott criterion are shown as black lines for the stripping radius as well as for the remaining gas mass. The gray line in 
the mass plots displays the predicted remaining gas mass if only continuous stripping is taken into account. In cases where the galaxy 
was not stripped completely, we did a linear fit for the bound gas mass for the final Gyr of the simulations as well as for the analytical 
estimate of continuous stripping. The resulting mass loss rates are denoted colour-coded in the plot. 
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between the minimum and maximum radius lines in Fig. [71 
A further agreement with the constant wind simulations is 
the fact that for the highest ram pressure peak, the galaxy 
is stripped completely even if it is moving near-edge-on (run 
Cl-SM-FST-E). 

The slightly premature complete stripping of the galax- 
ies in runs Cl-SM-FST-MF and C3-SM-FST-MF is caused 
by the fact that the remaining gas disc is affected by numer- 
ical diffusion due to the relatively low resolution of 0.5 kpc. 
Details are discussed in App.|A] All other results concerning 
ram pressure pushing are not affected by numerical diffusion 
or resolution. 

For the remaining gas mass, the simple analytical esti- 
mate and the numerical results agree only roughly. Gener- 
ally, in the simulations the galaxies lose their gas more slowly 
than predicted, because the analytical estimate does not in- 
clude the delay in mass loss explained above (Sect. 13.1)) . 
Moreover, it neglects effects of inclination and continuous 
stripping. 

3.3.2 Continuous stripping 

The gray line in the mass plots in Fig. [3 displays the analyt- 
ical estimate for the remaining gas mass if purely continuous 
stripping (see Eq. [4} is taken into account. Before pericen- 
tre passage, this is most likely to apply to galaxies moving 
edge-on. Indeed, we find a good agreement in mass loss rate 
for the first Gyr of runs Cl-LG-FST-EF and Cl-SM-FST-E, 
where this condition is fulfilled. For C2-SM-SLW-EMF, the 
ram pressure is very small during the first Gyr and thus the 
delay in mass loss is long. Therefore, the mass loss rate is 
overpredicted here. 

Pure continuous stripping is also expected to apply to 
the time well after pericentre passage. Thus, we have cal- 
culated averaged mass loss rates for the final Gyr of the 
simulations. We have also calculated the averaged mass loss 
rate for this time interval from the analytical prediction by 
Eq. 2] The values are given in Fig. [7] The mass loss rates 
at this late stage of the simulations have to be interpreted 
as upper limits, because due to the relatively low resolu- 
tion of 0.5 kpc the remaining gas disc is already significantly 
affected by numerical diffusion (see Sect. [A] for a detailed 
discussion) . 



4 DISCUSSION AND SUMMARY 

We have presented 3D hydrodynamical simulations of RPS 
of a disc galaxy orbiting in a galaxy cluster. We have focused 
on the evolution of the radius of the remaining gas disc and 
the mass of gas bound to the galaxy. In particular, we com- 
pared the numerical results to a time-dependent version of 
the classical GG72 estimate. 

For galaxies not moving near edge-on, we find that the 
stripping radius is predicted well by the estimate based on 
the GG72 criterion, unless the ram pressure increases faster 
than the stripping timescale. Then the remaining gas disc is 
somewhat larger than predicted. However, the evolution of 
the bound gas mass differs between simulations and the ana- 
lytical estimate: In the simulations, the galaxy loses the gas 
more slowly than predicted. Deviations of the analytical es- 
timate from the simulation results are caused by effects that 



are not taken into account by the estimate: the inclination of 
the galaxy, the time delay of true gas loss, and the combined 
effect of ram pressure pushing and continuous stripping. Al- 
though our simulations are an important step towards re- 
alistic models of RPS, they also neglect some aspects that 
have to be studied in the future, e.g. the multiphase physics 
of the ISM and transport processes and ambient motions in 
the ICM. 

The differences between analytical estimate and simula- 
tions in the evolution of the remaining gas mass will also be 
reflected in the distribution of stripped ISM throughout the 
cluster. If galaxies followed the prediction of the estimate 
based on the GG72 criterion, they would not lose gas af- 
ter pericentre passage. In the simulations, however, they do. 
Thus, the distribution of the stripped ISM along the orbit 
is broader than expected. 
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APPENDIX A: INFLUENCE OF RESOLUTION 

In order to check the influence of the resolution on our re- 
sults, we repeated two runs with better resolution. Run Cl- 
LG-FST-MF was repeated such that the resolution was im- 
proved by a factor of 2 everywhere in the grid compared to 
the description in Sect. ED Run C2-SM-SLW-MFM was re- 
peated such that the inner disc region, i.e. a cylinder around 
the galactic centre with radius 10 kpc and thickness 6 kpc, 
was always refined to 0.25 kpc. 

In grid-codes, advection of mass and momentum into 
non-axial directions is always accompanied by a numerical 
viscosity due to numerical diffusion. In the rotating gas disc, 
the gas has to move on circular orbits in a Cartesian grid. 
This geometrical mismatch is strongest in the inner part of 
the galaxy. Below a certain radius, the resolution of the ro- 
tational motion becomes insufficient and the numerical dif- 
fusion introduces radial velocities, which lead to a decrease 
of the central density. Fig. IA1I shows an example for this 
effect. As expected, the decrease in central density is much 
weaker if if higher resolution is used. Interestingly, it is also 
less pronounced for galaxies moving near edge-on. 

For the evolution of the galaxy's gas disc mass and ra- 
dius, two consequences arise from this effect: 

Firstly, in cases where the galaxy moves with medium 
to face-on inclination, not only the central density, but also 
the central projected surface density, E, decreases, though 
much milder. The decrease in surface density makes the gas 
disc more vulnerable to ram pressure pushing (compare to 
Eq. [2}. For quite strong ram pressures as in run C2-SM- 
SLW-MFM, this could cause the galaxy centre to be stripped 
prematurely. This indeed happened in the low resolution 
case of run C2-SM-SLW-MFM. In the high resolution case, 
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Figure Al. Radial profiles of the density, p, projected surface 
density, X, and total velocity, vtot, in the galactic plane at t = 
1300 Myr for different runs along the orbit C2-SM-SLW (see leg- 
end). One dot for each cell in the galactic plane is plotted. Black 
dots show the initial profile. Red and blue dots are for runs with 
inclination MFM, but different resolution (red for Ax = 0.25 kpc, 
denoted "HR" in the legend; blue for Ax = 0.5 kpc, denoted 
"LR"). Green dots are for inclination EMF and also resolution 
Ax = 0.5 kpc. 



the gas loss behaves as expected. Fig. IA2I shows the evo- 
lution of mass and radius for the two different resolutions. 
Up to pericentre passage, the two runs are nearly identical. 
The only other "critical" cases in the sense described above 
are the runs Cl-SM-FST-MF and C3-SM-FST-MF. For both 
runs, the gas disc radius and mass evolution agrees well with 
the analytical estimate up to shortly before pericentre pas- 
sage. In both runs, the galaxy is stripped completely shortly 
before pericentre passage, which is slightly faster than pre- 
dicted by the analytical estimate. We conclude that this sud- 
den gas loss shortly before pericentre passage is due to the 
numerical effect described above. However, we also conclude 
that the galaxies would indeed be stripped completely as the 
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Figure A2. Same as Fig. but comparing the high resolution 
(HR) and low resolution (LR) run for case C2-SM-SLW-MFM. 
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Figure A3. Same as Fig. IA2I but comparing the high resolution 
(HR) and low resolution (LR) run for case Cl-LG-FST-EF. 

one in run Cl-SM-FST-E is stripped completely but is not 
expected to suffer strongly from the numerical viscosity, as 
shown in Fig. IA1I 

The second consequence concerns continuous stripping. 
Figure IA3I compares the evolution of the mass and radius 
for run Cl-LG-FST-EF for two different resolutions. Up to 
~ 1500 Myr, both results agree very well. After that, the low 
resolution run yields a higher mass loss rate than the high 
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Figure A4. Colour-coded gas density in the orbital plane for run Cl-LG-FST-EF, at time t = 2100 Myr. The left panel is for high 
resolution (HR), the right panel for low resolution (LR). The black line marks the galaxy's orbit. 



resolution run. This is due to the fact that the influence of 
the numerical diffusion on the shape of the remaining gas 
disc is no longer negligible, as can be seen in 

Fig. El 

Concerning the evolution of the mass and radius of the 
remaining gas disc, we conclude that the results presented 
in Fig. [7] are reliable at least up to pericentre passage. After 
pericentre passage, the mass loss rates due to continuous 
stripping derived by our simulations, are upper limits. 
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